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(N ! Abstract 

o 

Finite field spreading scheme is proposed for a synchronous multiple-access channel with 
Gaussian noise and equal-power users. For each user, s information bits are spread jointly into 
a length-si vector by L multiplications on GF(2 ! ). Thus, each information bit is dispersed into sL 
transmitted symbols, and the finite field despreading (FF-DES) of each bit can take advantage of sL 
independent receiving observations. To show the performance gain of joint spreading quantitatively, 
an extrinsic information transfer (EXIT) function analysis of the FF-DES is given. It shows that the 
asymptotic slope of this EXIT function increases as s increases and is in fact the absolute slope of 
the bit error rate (BER) curve at the low BER region. This means that by increasing the length s 
of information bits for joint spreading, a larger absolute slope of the BER curve is achieved. For 
s, L > 2, the BER curve of the finite field spreading has a larger absolute slope than that of the 
single-user transmission with BPSK modulation. 

(N 
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I. Introduction 

In a Abuser multiple-access channel (MAC), each user regards the signals of other users 
as interference. When the number of users K is large, each user has a very low signal- 
to-interference-and-noise ratio (SINR) [1]. For this reason, spreading is usually employed 
as an SINR amplifier for each user, such as the conventional code-division multiple-access 
(CDMA) [2] and interleave-division multiple-access (IDMA) systems H-©. 

In the CDMA and IDMA systems, each information bit is spread independently. Specif- 
ically, in the CDMA system each information bit is spread by multiplying a binary vector 
into a length-L bit-vector, which is then interleaved by a bit-level interleaver. Since each 
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information bit has L independent receiving observations, the despreading (DES) will output 
an ameliorated signal with a sufficient large SINR for a further outer decoding (if there is 
a channel code as an outer code). In the IDMA system, a chip-level interleaving is jointly 
performed on multiple bit-vectors instead of the bit-level interleaving in the CDMA system 
HI. As the number of bit- vectors for interleaving is large, the IDMA in fact becomes a 
multi-user sparse-graph code that is appropriate for iterative decoding [6J. For this reason, 
the IDMA under an iterative chip-by-chip decoding provides a lower bit error rate (BER) 
than the conventional CDMA Q I0. Simulations in [3] showed that at the low BER region, 
the BER of the uncoded IDMA system can converge to that of single-user transmission 
with BPSK modulation. The BER of single-user transmission with BPSK modulation is the 
performance upper bound for both of the CDMA and IDMA systems, since in both systems 
the independent spreading of each information bit implies that each user employs a repetition 
code. 

In this paper, we propose a finite field spreading scheme for a synchronous MAC with 
Gaussian noise and equal-power users. For each user, every s information bits are spread 
jointly into a length-,?/, vector by L multiplications on GF(2*). A chip-level interleaving is 
then performed to generate the transmitted vector to the MAC. At the receiver, a multi-user 
iterative decoding is performed on a single factor graph to recover the information vector 
of each user. In our scheme, due to the joint spreading, each information bit is dispersed 
into sL transmitted symbols, and the finite field despreading (FF-DES) of each bit can take 
advantage of sL independent receiving observations. To show the performance gain of joint 
spreading quantitatively, we analyze the extrinsic information transfer (EXIT) function of the 
FF-DES. We show that the asymptotic slope of this EXIT function increases as s increases 
and is in fact the absolute slope of the BER curve at the low BER region. This means that 
by increasing the length s of information bits for joint spreading, a larger absolute slope of 
the BER curve is achieved. For s, L > 2, the BER curve of the finite field spreading has a 
larger absolute slope than that of the single-user transmission with BPSK modulation. 

II. Transmitter Model 

Figure \T\ illustrates a block diagram of the transmitter of ^-user finite field spreading 
multiple-access system. The length- .s^/V information vector u {k) = (u\ \u { 2 \ ...,u ( s ^),u\ ) e 
{+1, -1} = X, 1 < k < K, of user k is first mapped into a length-JV vector {k) = (J3f,pf, ...,$/) 
over GF(2 i ), i.e., every s information bits are mapped into a field element by mapping 
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Fig. 1. Transmitter of A'-user finite field spreading multiple-access system. 



r : X s — > GF(2 S ). Here GF(2 S ) = {0, 1, a..., a } is a finite field with a primitive element a, 



o<) 



?« ; c 



and T can be an arbitrary bijection from X s to GF(2 i ). Each element J3) ,1 < j < N, in fi ' is 



j 



spread by length-L spreading vector s® over GF(2") into J3 . s {k) = (J3) s) ,B\ si , ...,&, si ), 

J J * J *■ J *-* 

where the multiplication is on GF(2 V ). Here s (k) can be an arbitrary vector over GF(20 



with si' =£ 0, 1 = 1,2, ...,L. The output field vector after multiplication, denoted as y {k) = 



..(*) »,(*) 



.(*) 



,(*) Jk) 



Jk) 



(y\ ,y 2 , ...,y NL ), is demapped into binary vector c {) = (c\ ,c\ , ,.,,c4 I ), i.e., each field 
element is demapped into s bits by demapping T" 1 : GF(2 i ) — > X s , where T" 1 is the inverse 
transform of T. Vector c w , referred to as chip vector, is interleaved by a length- sNL chip- 
level interleaver n^ k) and is multiplied by amplitude J^- to generate the transmitted vector 
^f x (k) = (■ypjfxf, Jf xf,..., ^£xf NL \xf ] 6 X, to the Gaussian MAC. Here E h is the 
energy per information bit, and -£ is the symbol energy per transmission due to code rate £ 
of each user. The chip-level interleaving should be different for each user as in the IDMA 
system (3-10. 

The receiver receives a superimposed signal vector y = (yuyi, ...,y S Ni) with 



k=\ ' 



(1) 



where z t is a zero-mean Gaussian variable with one-side power spectral density Nq. Here, we 
assume that the symbols from the K users are synchronous. An iterative multi-user decoding 
is performed to recover information vectors u ik \k = 1, ...,K. 

To simplify the description, the mapping, multiplication of spreading vector, and demapping 
in the dotted box in Fig. Q] is referred to as finite field spreading. It should be noted that 
by the finite field spreading over GF(2 V ), the spreading of s information bits is performed 
jointly, and each information bit is dispersed into sL transmitted symbols. Thus, the FF- 
DES of each information bit can take advantage of sL independent receiving observations. 
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Fig. 2. Factor graph of the finite field spreading and the MAC for user k. 

Conventional spreading schemes, such as the CDMA and IDMA where the information bit 
is spread independently, are special cases of the finite field spreading with s = 1 . 

III. Iterative Decoding on Factor Graph 

Due to the chip-level interleaving, as length N is large, the ^-user finite field spreading 
multiple-access can be regarded as a Abuser sparse-graph code which can be decoded itera- 
tively on a single factor graph. In this section, we give an iterative decoding algorithm of the 
finite field spreading multiple-access on a single factor graph. In Section IIII-Al we represent 
the Abuser finite field spreading and the MAC by a single factor graph. In Section IIII-Bl we 
give the iterative decoding algorithm. 

A. Factor Graph 

Since each user has the same encoding process, the factor graph of each user is the same. 
We only illustrate the factor graph of user k in Fig. [21 There are three kinds of nodes in the 
graph: variable, mapping, and sum nodes. Let U, M a , V, Mb, and Y be five node sets that 
include the nodes below letters "[/," "M a ," "V," "M b ," and "F' in Fig. [2J The variable node 
in U corresponds to an information bit. The variable node in V corresponds to a symbol in 
GF(2"). The mapping node in M a or My denotes a mapping relation between a GF(2") symbol 
and s bits. The sum node in Y, associated with a received symbol, denotes a superposition 
of the transmitted symbols from K users. Note that the sum node connects the remaining 
K - 1 user's factor graphs. Edges connecting variable nodes in V and mapping nodes in Mb 
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are labeled by elements of spreading vector s (k) . During encoding and decoding, the message 
should multiply the labeled element or its inverse when passing across one of these edges. 

B. Iterative Decoding 

The iterative decoding is performed on the factor graph in Fig. [2] and is accomplished 
by efficient local decoding at all the nodes and interactions. The input and output of local 
decoding at each node is denoted by a log-likelihood ratio (LLR). Since the decoding of each 
user is performed in parallel and is the same, we only give decoding algorithm for user k. 

A single decoding iteration includes the local decoding at the sum node in Y, deinterleaving, 
and the FF-DES in Fig. [2] For the local decoding at the sum node in Y, we employ the low- 
complexity elementary signal estimation (ESE) algorithm given in [|3]|-[[5]|. For the FF-DES, 
we give a maximum a posteriori probability (MAP) decoding algorithm. After a number of 
iterations, a hard decision is performed to recover the information vector. 

1 ) Sum Node Decoding (ESE): Let L a (xf), 1 < i < K,i # k, denote a priori LLR of xf 
in (OQ) to the t-th, 1 < t < sNL, sum node in Y. By regarding J^f =[ i+k -J^-v as a Gaussian 
variable, this sum node performs a local decoding and outputs an extrinsic LLR of x\ } as 
0-0 

r(^)-iog Pr( ^ ) = +11 ^ 

L{Xt) l08 Pr(^ = -Hv,) 

^tanh(^)) 

(2) 



2jf(y f - Jf2f=i^tanh(^) 



f Z f =w (l-(tanh(^)p) + f 

which will be deinterleaved as a priori LLR for the FF-DES. 

2) FF-DES: After deinterleaving, we obtain a priori LLR L a (c t ), 1 < t < sNL, for each 
chip of user k. The FF-DES will calculate an extrinsic LLR L e (c t ) based on the priori LLRs 
L a (c\ '), 1 < t < sNL, by the local decoding at the nodes in M h — > V — > M h , subsequently. 

Since the decoding operation for each chip is the same, to simplify our description, we 
only introduce the calculation of extrinsic LLR L e (cL_ 1)s+;i ), \ < I < L, \ <n < s, performed 
on the subgraph associated with the first variable node in V of Fig. [21 

Let (L a (c ( (j l l)s+l ), ...,L a (c i h ) )) denote the s chip-LLR inputs to the ?'-th, 1 < i < L, mapping 
node in Mb and y\ = CV( r _ 1)s+1 , ..., y' is ) denote the associated received symbols iy' = (y[, ■■■,y' sNL ) 
is the deinterleaved form of y). This mapping node performs a local decoding to transform 
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a,AkK _ f T a,Ak\ r «/„ftK r« ^,(k),, , (k) 



chip LLRs to field symbol LLR vector L\yf) = (L a (y) K >), L\(yf), ...,L a aV _ 2 {yf)) of yf as 

^(y- = log - (t) - - 



= log 



ll m =i P r ( c (;-i).v+ m - r m (<fy\y[i_iy s+m ) 

n m =i Pr ( c (( -i), +m = r~ (Q)\y[i-i )s+m ) 

„_ Pr(c*, w =r m '(/i)|y;. IW ) 
ll«i=l nrfj) 



. ( C (i-l)s+m l y (i-l)s+m' 

° 8 n .v Pr ^-(f-i), +m = r » J 1 (Q)l^,-- 1) ., +m ) 

llm=l jw (*) = _ 1 , / -, 



-If 



(i+r„'(« 



,.<*) 



n.V / Pr ( C (f-l).;+m~ +1 lyQ-l).t+m) 

m=1 I Pr(c W =-llv' ) 

= log 



n.V I ^ C (i-l).i+m + ^'(/-Q.t+m-* 

m=1 I Pr(c W =-llv' 1 

\ y (i-\)s+m »{i-l)s+m> 

2 



^'W-^<°W > > , = ,i,..V- 2 



wi=l 



where T^{A) takes the m-th bit of the demapped vector T~\X). Here we use a posteriori 
probability of zero element as the denominator in a field symbol LLR. 

Let yL = (y[, ...jj_jjj +1 , ... ,y' L ) denote the joint vector of y' { ,i = 1, ...,L,i ± I. Similarly, 
we have denotation sj = (s\,...,s e _ v s i l 1 ,...,s L ) ). Based on the L— 1 field symbol LLR 
vectors L a (y) ), 1 < i < L,i ± t, the first variable node in V performs a local decoding to 
calculate the extrinsic LLR vector L\yf) = {L e {yf),L\{yf), ...,Z/ 2J _ 2 (yf )) of yf as 

Te ,(kK , Pr ^f : jyjg> 

Pr0f> = A(sf y^') 

= lo g ~^w 



Pr(^ = 0^) 

PrO 
log 



1 "VtfJ 



, n^Pr(yf = A(g)^^ 

= log 

nii^prcrf^oij;) 

L 

= Z W^' A = 0,l,...,or 2s - 2 (4) 

where 0l_i = (0, ...,0) is a length-(L - 1) zero vector and (4 V 1 is the inverse of s ( { } in 
GF(2 S ). 
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The ^-th, 1 < I < L, mapping node in Mb transforms the extrinsic symbol LLR vector 
U{yf ) ) to extrinsic chip LLRs of c® l)s+n as 



= log 



{l-\)s+n L \J"(t' 

Z J/ leGF(2'),r,7 1 (/i)=+l ^ T Wf = ^LV|f) 



LAeGF(2*) 2 Pr( 7 ?U)\y> e ) 
2^eGF(2*) 2 p r(y W =0[y j f ) 

, Z^gf(2-)(i +r; 1 (^))^ (r * )) 

= log — , n=l,...,s. (5) 

ZieGF^l-r,: 1 ^))^^ 

These extrinsic chip LLRs will be interleaved to update the priori LLRs of the sum node 
decoding in ©. 

3) Hard Decision: The hard decision is performed on the output from the mapping node 
in M a to the variable node in U. Using a similar principle as in ©, the first variable node in 
V calculates a total posteriori LLR vector L(fif ) ) = (L (J3\ k) ),L 1 (J3[ k) ), ...,L a v-2(fif)) of ySf as 

L 

The first mapping node in M a transforms this field symbol LLR vector to s bit LLRs using 

a similar principle as in © 

Zi e GF(2n(l+r; 1 (i))e L ' 1 < )) 
L(u n ) = log —-J-, n = 1, ..., s (7) 

for bit decision. Note that the mapping node in M a does not provide any message to the 
variable node in V during the iterative decoding. Here we give a bit decision algorithm, 
in which the hard decision is performed on the LLRs of information bits in ©. The hard 
decision can also be performed on the field symbol LLR in © to obtain an estimation of 

IV. Analysis of EXIT Function of FF-DES 

To show the performance gain of joint spreading quantitatively, we analyze the EXIT 
function of the FF-DES. In Section HV-Al we give an approximate EXIT function and show 
that this EXIT function asymptotically approaches a line. In Section IIV-BI we derive the 
asymptotic slope of the EXIT function. 
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A. Approximate EXIT Function 

The EXIT function describes a relation between the priori input and the extrinsic output of 
a decoding. Both the input and the output of the decoding are measured by mutual information 
or LLR mean value based on the Gaussian approximation [PTII- lfTOll . In this section, we give 
an approximate EXIT function of the FF-DES using the measure of LLR mean value. 

The EXIT function of FF-DES describes the relation between the mean value of a priori 
chip LLR and that of the extrinsic chip LLR. Generally, this EXIT function is determined by 
the realizations of a specific chip, mapping Y, and spreading vector s (k \ Here we consider 
random chip, random mapping Y, and random spreading vector s {k \ all of which are uniformly 
generated from all their possible realizations. We derive an expected EXIT function that is 
an average for all the possible chip, mapping, and spreading vector realizations. 

Since as stated in Section Unl the decoding operation for each chip of each user is the same, 
we omit superscript (k) to consider the chip %_!),,+„, 1 < £ < L, 1 < n < s, in our analysis. The 
EXIT function is the function between E[c^-\) S+n L e (c^-i) S+n )] and E[c(i-y )s+m L a (c(i-\) S+m )], i ± £, 
1 < m < s, where E[-] takes the expectation of a random variable. 

Combining ©, ©, and © in the FF-DES, we write L e (c {t - l)s+n ) as a function of L a (c {i - l)s+m ), 
i # £, 1 < m < s, 



L e (c {e - l)s+n ) = log 



^ „ „,,,„ yL y* r-,'(i( Jf )-'. ! ,)-r,-'(0) 

Z^GF^O-rVW)^ 1 - 2 "- 5 L (f< '- 1)!+ '" ) 



08 Z,i e GF(2o(i-r„- 1 (^))^ ( 

where p(A) = \ Yl l>m Y, s m= ^- m \A{s c r l sdL a {c {i ., )s+m ). 

By the Gaussian approximation (HI, C(,_i) iS . +m L"(c(,_i) i+m ) ~ N(m a ,2m a ), 1 < i < L, 1 < 
m < s, i.i.d., and c i (-Y )s+n L e {c i e-i) S + n ) ~ N{m e ,2m e ), where N(ju,o- 2 ) denotes the Gaussian 
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distribution with mean /i and variance a 2 . Based on dU), the EXIT function becomes 

m e = E[c(f-\) S+n L (c(f-i) S+n )] 

, Z. e GF(2 J )(i+r n 1 W)^ ) 1 
= jE[C( '- 1)Mlog Z^2,(i-^W)^ 1 

(l + c (f . 1)A+n r,; 1 (y f ))e p(yf) + Z^gf^W 1 + c (e -, )s+n T-\A))e^ 
og (l- c {e _ l)s+n T-Kyt))e<>w + Z^gf(2o,^(1- c^D^r^w)^^ 

2ieGF(2'),i# K (l - C(f-\)s+rF n l (X))ePW 

pn , ,, A pn 1 + ^. e GF(2^u^(i+%- 1) , +n r,; 1 (^))^ ) -^ ) i 

= E[2p(y e )] + E[\og — : ] 

i 2A6GF(20 ^(i_ C( ,_ 1>+n r n -i(^))^> + ^) 

= j(L - \)m a - E[\og(Yj e p(ye)+p{A) )] + £[log(l + J] g-WrdnPW)))]. (10) 

AeA- AeA+ 

Eq. © is from the fact Y~ l {j{) = c^-i)s+n since y c is the correct symbol. In (flOT ). the first 
term is from p(y^) = ± £f =1>! ^ ZLi C(.'-i)«+».^(c(J-i)i+«) due to r^fofo) -1 ^ = c ( ,-_i )5+m . 
In the second and third terms, A~ = {A\A 6 GF(2'),I^ 1 U) = -%-!),+„} and i + = {A\A e 
GF(2 S ),A * y c ,T~\A) = %_i )v+n }. It holds that \A~\ = 2 S ~\ \A + \ = 2 s " 1 - 1, where | • | takes the 
cardinality of a set. 

Now we give an approximation for the EXIT function in (flOl) . For A e A~ , we have 

1 L s 

p(ye)+p(X) = ~ Jj ^ i (c (i ^ )s+m +r~ m \A(s £ )~ l s i ))L a (c (i - 1)s+m ) 



2 

i=l,lW»!=l 



(1 +c {i -i )s+m r m \A(s e ) l Sj)) 



2 



L s n „ r-i/ hi c \-\* 



- ^, /_j ~ C((-_i) J+OT Z^ lC(j_i)j,.,. m J. U-U 

Similarly, for /I e A + , we have 

P(W-PW= 2j Zj 7 c (i _i), +m L (c ( ,_i), +m ). (12) 

Vector (C(,_i )s+ ir[ 1 (/l(5^)" 1 5,), ...jC^rj 1 ^^) -1 *,)) in (fTTI) and (PT2l) is a bit-wise correlation 
between the binary vectors for the correct symbol and an error symbol. Due to random chip 
C(j-i) S+m , random mapping Y, and random spreading element *,, this bit-wise correlation vector 
is approximately uniformly distributed on the set that includes all the binary vector in X s 
except all-one vector l s = (1, ..., 1). Let set Q,~ include all the binary-{0, 1} vectors of length 
s except l v , and f2+ include all the binary-{0, 1} vectors of length s except S . Let ti(Q.) 
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to 

denote the uniform distribution on set Q.. We have the following approximation of the EXIT 
function CGI 

2 s - 1 T 2 r -'-l T 

m e ^5(L-l)m a - J B[log(^e^l 1 ^) ]+jB[log(1+ £ e-^' r ^)]^ w (13) 

7=1 7=1 

where rp,j = 1, ...,2 s ~ 1 ,i = 1,...,L-1, are i.i.d. with rj >i ~ < U(Q.~). Vectors r' ., j=\, ...,2 V ~ 1 -1, 
z = 1,...,L- 1, are i.i.d. with r'.. ~ 1/(Q+). Elements /?,>,/ = l,...,L,m = l,...,*, in ht,i = 
1, ...,L, are all i.i.d. with h im ~ N(m a ,2m a ). 

To see the accuracy of the approximation in (|T~3T >. Figs. |3] and |4] illustrate the EXIT function 
in (flOl) and y>(m a ) in ([131) for L = 8, 16 and 5 = 1,2,4,6 by Monte Carlo simulations. We 
see that <p(m a ) is in fact a tight upper bound of the EXIT function. The ratio between (p(m a ) 
and the EXIT function is less than 1.06 for m a < 10 in both figures. Based on this fact, in 
the rest analysis of this paper, we employ approximate EXIT function (p(m a ) of the FF-DES. 

Moreover, we have the following observation on the approximate EXIT function. There 
exists a cutoff point m usually less than 1. For m a > mo, every curve of the approximate 
EXIT function in Figs. [3] and |4] asymptotically approaches a line. The asymptotic slope of 
the approximate EXIT function increases as s increases. This means that by joint spreading 
for s information bits, an SINR gain is obtained with respect to that of s = 1. This SINR 
gain will be enhanced as length s increases. Furthermore, by comparing Fig. |3]and Fig. @]we 
can see that given s, the asymptotic slope of the approximate EXIT function for L = 16 is 
larger than that for L = 8. This is due to that as spreading length L increases, the advantage 
of joint spreading is enlarged. 

Remark 1: Although we have focused on the asymptotic slope of the EXIT function of 
the FF-DES when m a is large above, the EXIT function at m a near is also noteworthy. 
If the number of users K is large, due to a large multi-user interference, at the beginning 
of decoding iteration the extrinsic output of the ESE is always very small. This requires 
the EXIT function of the FF-DES to be sufficiently large at m a near 0, or else the iterative 
decoding would fail at the beginning of iteration. As revealed in IfTTTl and lfT2l . conventional 
spreading scheme is widely used in the MAC since it can provide a large EXIT function at 
m a near 0. We see that in Figs. [3] and @] the finite field spreading can provide a similar EXIT 
function as that of the conventional spreading scheme (s = 1) at m a < m , and can also work 
in the environment of large multi-user interference. 
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Fig. 3. EXIT functions (solid lines) and their approximations of <p(m a ) (dashed lines) of the FF-DES for L = 8 and 
s= 1,2,4,6. 
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Fig. 4. EXIT functions (solid lines) and their approximations of <p(m„) (dashed lines) of the FF-DES for L = 16 and 
s= 1,2,4,6. 



B. Asymptotic Slope 

Section HV- Al has shown that when m a is larger than cutoff point m , the approximate EXIT 
function of FF-DES asymptotically approaches a line. In this section, we analyze this slope 
by deriving the asymptotic differential of (p(m a ) in ([131) . 

Theorem 1: Let g(s,L) be a function of s and L. Let m > be a constant. Suppose that 
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12 



<p(m a ) = <p(tno) + g(s, L)(m a - m ) holds for m a > mo, we have 



(,v-l)(L-l) 



( L-\ i v \ 



^ l ) =l - i+ (2 ,_ 1)( ^- z z np 

V ' k=\ \(n u ...,n,- ] )e0(k-\) i=l v l > 



2-1 



(14) 



where #(fc-l) - {m, ...,nz,_i| 2^ ft; < fc-1,0 < rc,- < 5-l,n,- e Z} with integer set Z. □ 

To prove Theorem [1] we first prove the following lemma. 
Lemma 1: Let a > 1 be a constant. Let r} = |/}i, ...,T ;m ],T ; -,- 6 {—1,0, 1 }, j = 1, ...,«, be 

i.i.d. random vectors with w(r}) - 2/li O,' - ^» an( ^ W+ (T/) - 2Xi (r/.»'l > ^ f° r a ^ ./'• Let ^ 
be a length-m vector whose elements are i.i.d. with Gaussian distribution Nip,, 2//),// > 0. It 
holds that T 

r E[l0g(a + ^ ^ )] n 

lim = 0. (15) 



D 



Proof: We first have 
< £[log(a + V e r ^ )] < £[log(a + V /^ -"O)] 

_ — T 

< £[log((n + l)max{a,e irih ~ w{rM , ...,e ir " h ~ w( - r M})] 

n 

< log(w + 1) + E[log a + J^ \Fjh - w(?j)fi\] 



7=1 ;=i 

^T — T 



7=1 
-T 



log(a(n + l))+/»£E[|ri» -w(?i)/i| |w + (ri) = fc]Pr(w + (?i) = fc) 



fc=l 

in 



Ikfx 



\og(a(n + l))+2n V J —Fr(w + (r l ) = k) (16) 

f-^ V 7T 



A-l 




< log(a(n + l))+2n A / — 

— T 

where ([TBI is due to that under the condition w + (?i) = k, (F\h -w(T])ij) ~ N(0,2k/j) holds 
for each realization of ri. Thus, 



~T 



£[log(a + 2"-i ^ )] lo g(«(" + 1 ))+ 2 " V^ 
< lim < lim — = 0. 

The lemma is proved. 

Proof of Theorem\l} Let j max - {j*\j* = argmax ; {£^r 1 1 w{r^)} be the set of subscripts that 
maximize the Hamming weight of joint vector (r^i, ...,rj nL -i) in (TT3T) . Given an element of 
7* e Jmax, we have set / = {j\r p = r /;! -, i = 1, ..., L - 1 } C j max . 
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Since (p(m a ) = (f(mo) + g{s,L){m a - m ) holds as m a — » +00 by assumption, using (PT3T) we 
have 

, n ,. <p(m a ) - <p(m ) (f(m a ) 

g(s,L) = lim = lim 

m a -*+oo m a — m$ ma->+oo m a 

£[log(ZS ^' r ^')] ^tlogd + Z^" 1 e-^nti)] 
= s{L - 1) - lim + lim 



m n 



m„ 



= s(L - 1) - lim 

m a — »+oo 



^[log(2S^- r ^0] 



m„ 



(17) 



^(L - 1) - lim lim 

m a ->+oo m a m a — >+oo 

L-l 



£[log(|/| + ZS 1 J , / e-^ 1 ^--^^)] 



m„ 



(18) 



r'=l 
Cs-1)(L-1) L-l 

= s(L-l)- £ fcPr(^w(r A ,) = fc) 

fc=l i=\ 

(JT-1)(L-1) L-l 

= <l - 1) - 2 ^ pr <Z w{r ^ d - k) - Pr( Z w(r ^- ;) - * " 1}) 

k=\ i=\ i=l 

L-l (*-l)(i-l) L-l 

= s(L-l) - (s-l)(L-l)Pr(J]w(r r j) < ( s -l)(L-l)) + J] P^w^j) * k-1) 

r=i k=\ 

(s-l)(L-l) L-l L-l 

=l-i+ Yj v<Yj w{ri < i) - k - l —Yj w{r2si > i) - k - l) 



L-l 



k=\ (=1 



k=\ i=\ 

(s-\)(L-\) / L-l 



\ 



i=\ 
2 s - 1 



= L-\+ Yj PrC^^u)^*" 1 ) 

jfe=l V 1=1 

Cs-1)(L-1) / ,.-1 , N 

t z z n : 

fc=l Urn ,...,/!,._ ,)e0(fc-l) 1=1 v J ' 



(19) 



L-l + 



(2 s - iy L - 



L-\ / \\ 



(20) 

1 1. i n . 1 
V(ni,...,«L-i)e6»(i-l) 

Eq. (fTTT) is from Lemma \T\ Eq. (fTSl is due to that elements of /i, are i.i.d. with distribution 
N(m a , 2m a ) and LemmaQ] Eq. <Qj3> is due to that Ym=\ w(fjj), 7=1, •••, 2 s " 1 , are i.i.d.. Eq. © 
is due to that r 1( , / = 1, ...,L- 1, are i.i.d. with ry ~ 1l(Q.~). The theorem is proved. n 

For s = 1, g(l, L) = L- 1 is the slope of the EXIT function of the DES for the conventional 
spreading scheme of IDMA [3]|-[[5]|. The last term in (1201) is the slope gain from the joint 
spreading. 

V. Slope of BER 

In this section, we show that the asymptotic slope of approximate EXIT function analyzed 
in Section IIV-BI in fact is the absolute slope of BER curve at the low BER region. We 
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verify our analysis by BER Monte Carlo simulations for practical finite field multiple-access 
systems. 

Let 4>{m a , ^0 be the EXIT function of the ESE (Appendix lAl). The mean of extrinsic LLR 
converges to infinity as the number of iteration increases if and only if <p(x, jj*-) > ^ -1 (*) holds 
for x > 0-01511, where ^ _1 (0 is the inverse function of <p(-) in (fl3T) . Since <p(x, ||) < 4^ 
(Lemma [2] in Appendix lAl) and lim A .^ +00 ^ _1 (jc) =+oo, to converge to infinite mean of LLR, 
Y should approach infinity. This is due to that no channel code is employed for each user. 

Since for a given x, ip~ l {x) is a decreasing function of spreading length L, <p(x, jf-) > tp~ l (x) 
holds for x > as L is sufficiently large and j^ —> +oo. In this case, the mean of the extrinsic 
LLR converges to infinity as the number of iteration increases. Due to lim A ^ +co cp(x, jf-) = 4^- 
from Lemma |2l the asymptotic BER of hard decision at the information node in U is estimated 
as 0-001 

P e * Q{ Jlg(s,L+\)-^-) < <r ?( *'< (21) 

where we have used tp(x) « g{s,L)x as x is large, Q{x) < e~ x ~ 12 , and ~g(s,L) = g(s,L+ \)/L 
is referred to as a standard slope. Note that we have used slope g(s, L+ 1) instead of g(s,L) 
in (|2TI) since the hard decision is performed based on the total LLR in © other than the 
extrinsic LLR in ©. 

In fact as the BER approaches 0, the upper bound in (I2TI) is very tight. Eq. (I2TT) indicates 
that at the low BER region, the BER curve is approximate to a line with absolute slope 
g~(s,L). Since ~g(l,L) = 1 holds for an arbitrary L, at the low BER region, the BER curve of 
conventional spreading scheme (s = 1) always converges to that of (single-user) uncoded 
BPSK transmission. Unfortunately, this predicament can not be improved by increasing 
spreading length L. In our work, since ~g(s, L) > 1 holds for s,L>2, the BER curve of finite 
field spreading with s,L>2 has a larger absolute slope than that of single-user transmission 
with BPSK modulation, and the absolute slope of the BER curve increases as s increases. 
In addition, there is a special case of L = 1, i.e., code rate of each user is 1. It is easy to 
show that in this case, the BER can never approach for the number of users K > 2 due to 
code rate 1 of each user. For single-user (K = 1) finite field spreading multiple-access, since 
~g(s, 1)^1 for an arbitrary s, at the low BER region, the absolute slope of the BER curve 
approaches that of single-user transmission with BPSK modulation. 

We verify our analysis above, obtained based on the Gaussian approximation and the 
approximation of (PT3T ). by BER Monte Carlo simulations on practical systems. We consider 

February 26, 2013 DRAFT 



15 










— L 


= 8 


-1. 




-l-L = 16 


"H«=*-~ T" ^^V 


-3 
-4 
-5 
-6 

5 -8 
-9 
-10 
-11 
-12 
-13 
-14 
.is 




\ At \ H 











Al\ 


\\ single user \ 










M \ 


s = V, 


i = r 






- 6 Wi4 






! I I ! \ 







5 6 



7 8 9 10 



Fig. 5. BER curves of (K = 8)-user finite field spreading multiple-access systems with L 
information vector length is sN = 12000 and the number of decoding iteration is 50. 



8,16 and s =1,2,4,6. The 



(K = 8)-user finite field spreading multiple-access system with L = 8, 16 and s = 1,2,4,6. 
Fig. [5] illustrates BER curves obtained by Monte Carlo simulations. Since the EXIT function 
theory is based on the assumption of infinite code length [TTII- lfTDl . in our simulations we 
employ a large code length with the information vector length of sN = 12000. The number of 
decoding iteration is 50. Both mapping T(-) and the chip-level interleaving are random, and 
the BER curve illustrates an average BER for all the possible mapping and interleaving real- 
izations. We use horizontal coordinate jj 1 and vertical coordinate log P e (natural logarithm). 
We see that at the low BER region, all BER curves are approximate to lines. We compare 
the absolute slopes of BER curves at low BER region in Fig. [5] with standard slope ~g(s,L) 
obtained by Theorem \T\ in Table HI The difference is less than 0.1 for all the comparison 
pairs. The two curves for s = 1 overlap with each other and converge to the BER curve 
of single-user transmission with BPSK modulation. Curves for s > 2 have larger absolute 
slopes than that of single-user transmission with BPSK modulation. All these phenomenons 
coincide with our analysis. 
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TABLE I 
Comparison between absolute slopes of BER curves at low BER region in Fig.[5]and ~g(s, L) 



L 


8 


16 


s 


1 


2 


4 


6 


1 


2 


4 


6 


BER slope 


1.0625 


1.2262 


1.6364 


2.2352 


1.0028 


1.2373 


1.7269 


2.3543 


g(s,L) 


1 


1.2411 


1.7002 


2.2095 


1 


1.2675 


1.8240 


2.4493 



VI. Conclusions 

In this paper, we proposed a finite field spreading scheme for the MAC, in which the 
spreading for multiple bits is performed jointly by the finite field multiplication. Under the 
iterative decoding, an SINR gain is obtained during the joint spreading. The multi-user finite 
field spreading multiple-access achieves a lower BER than that of the single-user transmission 
with BPSK modulation under the iterative decoding. 

We considered the finite field spreading multiple-access system without channel coding 
for each user. If there is a channel code employed as an outer code for each user, our EXIT 
function analysis of the FF-DES is also available from the theory of concatenated code 

EMU. 

There are still a number of problems in connection with our work that seem to deserve 
further investigation. In our analysis, we considered an average EXIT function and asymptotic 
slope for all the mapping and interleaving. In fact, different mapping and interleaving have 
different EXIT functions that provide different BER performance. Design of spreading and 
interleaving to achieve better BER performance is an interesting future work. 

In the field field spreading of Fig. El every s information bits are spread into L length- 
s' vectors, each of which can be regarded as a codeword of a rate-1 code. These L rate-1 
codes constitute a finite field spreading code with rate j. The L multiplications on finite field 
can be regarded as L rate-1 encoding operations that determine the L rate-1 codes. These 
rate-1 encoding can also realized by other methods, or we can employ codes with rate less 
than 1 to achieve a better EXIT transfer performance. All these problems deserve further 
investigations. 
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Appendix A 

EXIT FUNCTION OF ESE 

Using CQ) and ©, the EXIT function of the ESE is derived as 

r 2 Vf^( y ,-VfZ^tanh(^ )) 

fZf^(l-(tanh(^)) 2 ) + f 



^[x^rCxf)] = £[— "—— — -] 



2 f + 2f ^ Z£ w (^ - tanh(^)) + 2 ^ X fh t 

~ E[ fZr=u rt (l-(tanh(^))^) + f ] 

4 
2 Zf =urt (l - (tanhC^^ 2 ))^ + L/| 

= g[ ,„r-, 1 ..J " *0*. TT) 



22£ 1 (1 -(tanh(^))2) + L/| 'AT, 



o 



where »,■,* = 1, ...,K- 1 are i.i.d. with % ~ A/(f , f ). Eq. <E3 is due to that x^ ~ WO 
is independent of x\ , i ± k, and z t - Note that <f>(m a , fr) is an average EXIT function of the 



ESE for all the possible transmitted symbols. 

No' 



We give an upper bound for 4>{m a , ^) 



Lemma 2: 



tLu Lu 

<f>(m a ,-^)<4- 



with equality if m a — > +oo. □ 

Proof: The upper bound is from the fact |tanh(x)| < 1. We prove the equality in the upper 
bound. Let j < p < 1 be a constant. We have 

lim 0(m a ,^) > lim —^-^L \' Ka) J — (23) 

Wo 



11111 \1S\I I Lst^ I £_ 11111 _, 

m fl - + oo N »^+oo 2 2f = i 1 (l-(tanh(m a -(m fl )P)) 2 ) + L/f 



4C1 wa_^-l 

> lim — — (24) 

— + °° 2 Zf = 7 J (1 - (tanh(m fl - (m fl )/>)) 2 ) + L/g 



£ 



6 



where (1231 is from monotone increasing of tanh(jc) for x > 0, and (1241 is from Chebyshev 
inequality. Using the squeeze theorem we obtain lim ma ^ +0o <p(m a , fp) = 4^-. The lemma is 
proved. 
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